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We demonstrate that massive fields, such as dark matter, can directly produce a cosmological 
evolution of the fundamental constants of Nature. We show that a scalar or pseudoscalar (axion-like) 
dark matter field (f, which forms a coherently oscillating classical field and interacts with Standard 
Model particles via quadratic couplings in (f, produces ‘slow’ cosmological evolution and oscillating 
variations of the fundamental constants. We derive limits on the quadratic interactions of (f with the 
photon, electron and light quarks from measurements of the primordial ‘*He abundance produced 
during Big Bang nucleosynthesis and recent atomic dysprosium spectroscopy measurements. These 
limits improve on existing constraints by up to 15 orders of magnitude. We also derive limits on 
the previously unconstrained linear and quadratic interactions of cf with the massive vector bosons 
from measurements of the primordial "^He abundance. 
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Introduction. — The idea that the fundamental con¬ 
stants of Nature might vary with time can be traced as 
far back as the large numbers hypothesis of Dirac, who 
hypothesised that the gravitational constant G might be 
proportional to the reciprocal of the age of the Uni¬ 
verse [T]. More contemporary dark energy-type theo¬ 
ries, which predict the cosmological evolution of the fun¬ 
damental constants, such as Brans-Dicke models, string 
dilaton models, chameleon models and Bekenstein mod¬ 
els, assume that the underlying fields, which give rise to 
this evolution, are either massless or nearly massless, see 
e.g. Refs. [2Hl6]. The evolution of the underlying dark 
energy-type field is determined by its couplings to mat¬ 
ter, including dark matter, see e.g. Ref. [7]. Another 
possible way of achieving a variation of the fundamental 
constants is via quantum effects induced by cosmological 
renormalisation group flow, see e.g. Refs. [nno]. 

In this letter, we demonstrate that a cosmological evo¬ 
lution of the fundamental constants can arise directly 
from a massive dark matter (DM) field, which is not nec¬ 
essarily unnaturally light. The possibility of exploring 
DM models in this particular context opens an exciting 
new avenue in the study of the cosmological evolution of 
the fundamental constants, since DM models are more 
amenable to theoretical and experimental investigation 
compared with their dark energy-type counterparts (see 
e.g. Ref. [51] and references therein). 

One of the leading candidates for DM is the axion, 
a pseudoscalar particle which was originally introduced 
in order to resolve the strong CP problem of Quantum 
Chromodynamics (QCD) (55][53] (see also [531 - 157] ). The 
axion forms a coherently oscillating classical field in the 
early Universe, which survives to the present day if the 
axion is sufficiently light and weakly interacting. These 
relic axions, which now reside predominantly in the ob¬ 
served galactic DM haloes, can be sought for through a 
variety of distinctive signatures, see e.g. Refs. [581136] . 
Apart from the QCD axion, one can also consider other 


spin-0 DM candidates, including scalar DM, the effects 
of which in astrophysical and cosmological contexts have 
been investigated extensively, see e.g. Refs. [37H17]. 

In the present letter, we consider a non-relativistic cold 
scalar or pseudoscalar (axion-like) DM field (j>, which is 
produced non-thermally (e.g. through vacuum decay gsi- 
IMl) and forms a coherently oscillating classical field, cj) = 
(fo cos{ujt), that oscillates with frequency uj ~ m^c^/h, 
where is the mass of the DM particle, c is the speed of 
light and h is the reduced Planck constant m- In partic¬ 
ular, although {(/)) =0, ((^^) = (^q/2 7^ 0 for such a field. 
The DM model we consider in our present work satisfies 
constraints on both interaction strength and mass from 
existing experiments, including fifth-force searches and 
supernova energy loss bounds, and also satisfies gravita¬ 
tional requirements (including the formation of observed 
galactic DM haloes), since the non-gravitational interac¬ 
tions we consider are very weak. 

Theory. — The field (j) can couple to the Standard 
Model (SM) fields via the following quadratic-in-()) inter¬ 
actions [51] : 
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where the first term represents the coupling of the scalar 
field to the SM fermion fields /, with uif the standard 
mass of the fermion and / = the second term rep¬ 

resents the coupling of the scalar field to the electro¬ 
magnetic field tensor F, and the third term represents 
the coupling of the scalar field to the SM massive vector 
bosons V, with My the standard mass of the boson. Un¬ 
like the corresponding linear-in-t/) interactions, here (f can 
represent either a scalar or pseudoscalar particle, since p 
is of even parity for both. 
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Comparing the terms in Eq. Q with the relevant terms 
in the SM Lagrangian: 

^SM =-2_^mfff - - -: ( 2 ) 
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we see that the SM particle masses are altered as follows: 
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In order to see the effect of the quadratic coupling of (j) to 
the electromagnetic field tensor, it is convenient to write 
the relevant terms in an alternate system of units: 


Hubble friction. Thus, the temporal evolution of and 
spatial variations in produce ‘slow’ space-time vari¬ 
ations in the fundamental constants, which can be con¬ 
strained from astrophysical phenomena, most notably 
Big Bang nucleosynthesis (BBN) and cosmic microwave 
background (CMB) measurements [58], while the oscil¬ 
lating component of produces oscillating variations in 
the fundamental constants, which can be sought for with 
high-precision laboratory measurements, such as atomic 
clock and laser interferometry experiments |591160j . 

The energy density of a non-relativistic oscillating DM 
field is given by p ~ and evolves according to 

the relation: 

Pdm = 1.3 X 10-6 [l + z(t)]"^, ( 6 ) 
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from which we deduce that the electromagnetic fine- 
structure constant a is altered as follows: 
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provided that changes in (/> are adiabatic. 

Supernova energy loss arguments constrain the inter¬ 
action parameters that appear in Eq. Q, for the scalar 
masses < Pme 30 MeV. Consideration of the pho¬ 
ton pair-annihilation channel 7 -|- 7 —> nucleon 

bremmstrahlung channel A^-|-A^ -^N + N + 4) + (j) yield 
the limits: A! ^ 3 x 10 ^ GeV and > 15 x 10 ^ GeV, 
respectively 155] . The limits on the interaction parame¬ 
ters in Eq. Q from fifth-force searches are weaker, since 
for the quadratic couplings in Eq. Q , a fifth-force is pro¬ 
duced in the leading order by the exchange of a pair of 
^-quanta between two fermions, which generates a less 
efficient V{r) ~ —Trij:/64Tr^(A'p ■ Xjr^ attractive poten¬ 
tial, instead of the usual Yukawa potential in the case of 
linear couplings. The strongest fifth-force limits are for 
the proton interaction parameter [55]: A^ > 2 x 10^ GeV, 
for the scalar masses < 10“^ eV [55] . 

A lower limit on the mass of generic spin-0 particles, 
which saturate the observed cold DM content, comes 
from the requirement that the de Broglie wavelength of 
these particles not exceed the halo size of the smallest 
dwarf galaxies (i? ~ 1 kpc): > 10“^^ eV (however, 

this limit is relaxed if spin -0 particles are not the dom¬ 
inant contributor to cold DM). This simple estimate is 
in fact in good agreement with more sophisticated astro- 
physical and cosmological considerations [i^MT] . 

When ^ H{t), where rurj, is the mass of the 
DM particle and H{t) ~ l/2t is the Hubble param¬ 
eter as a function of time in the early Universe [57], 
(j) = 4>o cos{m^t) is an oscillating field and so (jP' contains 
both the oscillating term, (j)QCOs{2m^t)/2, as well as the 
non-oscillating term, = (()q/ 2. When H{t), 

(/) is a non-oscillating constant field due to the effects of 


where z{t) is the redshift parameter and the present mean 
DM energy density is determined from WMAP measure¬ 
ments [5T] (for relativistic DM, the mean DM energy 
density evolves as ppM oc [1 -|- z{t)]‘^). The present-day 
cold DM energy density in our local galactic region is 
/^CDM ~ 0.4 GeV/cm^ [5T]. The energy density of a non¬ 
oscillating DM field is given by p ~ /2 and, due 

to Hubble friction, is approximately constant while the 
field remains non-oscillating: 

Pdm — 1.3 X 10-6 [1 -I- z(tm)]^ -V, (7) 

cm'^ 

where z{tm) is defined by H(tm) ~ In both cases, 
the largest effect of variation of the fundamental con¬ 
stants induced by (/>, therefore, occurs during the earliest 
times of the Universe. 

CMB Constraints. — Variations in a and me at 
the time of electron-proton recombination affect the ion¬ 
isation fraction and Thomson scattering cross-section, 
O’Thomson = Stt/S • changing the mean-free-path 

length of photons at recombination and leaving distinct 
signatures in the CMB angular power spectrum. Recom¬ 
bination occurs over a relatively short period of time, 
AtcMs/^CMB ^ 1 ■with 2 ;(tcMB) ~ 1100 , meaning that 
the energy density of DM is approximately constant 
during recombination. Analysis of WMAP measure¬ 
ments, which give the bounds [lAa/a)cMB ^ 0.01 and 
(Ame/me)cMB ^ 0.04 [55], hence immediately yield the 
following bounds on the quadratic interactions of the (os¬ 
cillating) scalar field (p with the photon and electron: 
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where we have made use of the fact that the mean DM en¬ 
ergy density at the time of recombination is much greater 
than the present-day local cold DM energy density, and 
assumed that scalar DM saturates the present-day DM 
content. These constraints are presented in Fig. [^ 

BBN Constraints. — Variations in the fundamen¬ 
tal constants and particle masses in the early Universe 
alter the primordial abundances of the light elements. 
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The most sensitive constraints come from consideration 
of the primordial abundance of "^He. There are two dis¬ 
tinct regions to consider — the first is when the field 
(j) is non-oscillating and approximately constant dur¬ 
ing BBN, which corresponds to the scalar mass range 
TO 0 10“^® eV, while the second is when the field (j> is 
oscillating for the entire duration of BBN, which corre¬ 
sponds to the scalar mass range 10“^® eV. 

We begin with the former case, when (j> is constant. To 
leading order, changes in the primordial ^He abundance 
are given by: 


The terms on the right-hand side of Eq. 0 can be 
expressed in terms of variations of various fundamental 
constants: 
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where the variations in all of the parameters are at the 
time of weak interaction freeze-out, and: 
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where (n/p)weak = e^'^np/^woak jg neutron-to-proton 
ratio at the time of weak interaction freeze-out {Qnp = 
rUn — nip = aaAqcD + (nid — niu) is the neutron- 
proton mass difference, with the present-day values 
(aQ:AQCD)o = —0.76 MeV, where Aqcd ~ 250 MeV 
is the QCD scale and a is a dimensionless constant, 
and {md — to„)o = 2.05 MeV [53], while Tweak = 
biM^^ {6w )/« 0.75 MeV |S7| is the 

weak interaction freeze-out temperature, 0w is the Wein¬ 
berg angle, Mpianck is the Planck mass and bi is a dimen¬ 
sionless constant), r„ is the neutron decay rate given by 

(l/r„ = Tn^ 880 s gU): 


r„ = 


b2a'^ml 


sin"‘(0w) [15 
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+ x\n{x + \/x"^ — 1) , ( 10 ) 


with X = Qnpjnie and 62 a dimensionless constant, and 
tBBN ~ 180 s is the end-time of BBN. 
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where we have made use of the relation cos( 0 'w) = 
Mw/Mz, with (Mw/Mz)o = 0.882 [61]. In order to 
estimate the term AteBN/^BBN; we note that teBN oc 
Afpianck/TBBjj is determined from the condition that the 
rate of expansion and rate of strong interaction processes 
involved in BBN are equal. On dimensional grounds, we 
write Tbbn = b3-^^QCD/^pii,nck ~ ^0 keV, where 63 and y 
are dimensionless constants. For 63 ~ 1, we find y « 0.2, 
which gives: 
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From the measured and predicted (within the SM) pri¬ 
mordial "‘He abundance, l^®’"^P('‘He) = 0.2477 ± 0.0029 
[51] and l)(''®°’'(^He) = 0.2486 ± 0.0002 [55||55], we find 
the following constraints on the quadratic interactions of 
(j) with the SM particles when m^ <C 10“‘® eV, using 
Eqs. 0, @ and @: 
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~ (-0.5 ± 1.7) X 10"^° eV'‘, (14) 


where we have made use of the relation [1 + z(tm )]/(1 + 2 :weak) = V^weak/^m and the fact that the mean DM energy 
density during BBN is much greater than the present-day local cold DM energy density, and assumed that scalar DM 
saturates the present-day DM content. Here the k'x = ±1 correspond to the relevant signs in the Lagrangian Q. 
These constraints are presented in Fig. Since the scalar field is non-oscillating during BBN when m^ ^ 10“^^® eV, 
we likewise also have the following constraints on the analogous linear interactions of (j) with the SM particles: 
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We now consider the case when >> 10 eV, for 


which (j) is oscillating during BBN. In this case, the only 
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required modification in relation ([^ is the replacement 


A(r„tBBN) 


The dominant contribution to the variation of this in¬ 
tegral comes from the variation of the integrand, which 
has the time dependence Ar„(t) oc and from the 

variation of tweak- Changes in the primordial ^He abun¬ 
dance due to changes in the neutron lifetime are thus 


where tweak 


1.1 


s. 


suppressed by the small factor 2tweak/tBBN ~ 10“^ when 
(j) is an oscillating field, compared with the case when 
(j) is non-oscillating. From the measured and predicted 
(within the SM) primordial ^He abundance, we hence 
find the following constraints on the quadratic interac¬ 
tions of (j) with the SM particles when S> 10“^® eV, 
using Eqs. (§, ( [ 11 ] ), ( [T^ and the modified version of 
Eq. 
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where we have made use of the fact that the mean DM 
energy density during BBN is much greater than the 
present-day local cold DM energy density, and assumed 
that scalar DM saturates the present-day DM content. 
These constraints are presented in Fig. 

Atomic Spectroscopy Constraints. — We can de¬ 
rive constraints on the quadratic interaction of a coher¬ 
ently oscillating scalar field 4> with the photon through 
the oscillating shifts it induces in atomic transition fre¬ 
quencies that are sensitive to variations in a. Using the 
recent atomic dysprosium spectroscopy data of Ref. 
which were used to place constraints on the analogous 
linear interaction of (j) with the photon, we derive con¬ 
straints on the quadratic interaction of (/> with the pho¬ 
ton, assuming that scalar DM saturates the present-day 
local cold DM content. These constraints are presented 
in Fig.[TJ 

Spatial Variations in — We can derive con¬ 

straints on spatial variations in from measurements 


of spatial variations in the abundance of primordial deu¬ 
terium, as determined from quasar absorption spectra: 
d{D) = (5.4±2.9) X 10“®/ Glyr [SH] (we use the notation 
d{X) to denote the fractional spatial gradient in parame¬ 
ter X). Without resorting to a full, time-dependent, nu¬ 
merical calculation, we assume that spatial variations in 
the deuterium abundance are due to spatial variations in 
at the time of weak interaction freeze-out (for other 
mechanisms of generating spatial fluctuations in the pri¬ 
mordial light elemental abundances, see Refs. EH HO]): 
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where Zweak ~ 3.2 x 10® and 2 :, 
scalar DM saturates the present-day DM content, and 
using Eqs. (Iq and (11), this gives for » 10“®® eV: 
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Likewise, using Eqs. § and ((IT|), we find for 10 ®® eV: 
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where we have made use of the relation [1 -|- z{tm)]/{y + Zweak) = y/tweak/tm- 


Conclusions. — We have demonstrated that mas- sive fields, such as DM, can directly produce a cosmo- 
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logical evolution of the fundamental constants. We have 
shown that a (pseudo) scalar DM field (j), which forms 
a coherently oscillating classical field and interacts with 
SM particles via quadratic couplings in </>, produces ‘slow’ 
cosmological evolution and oscillating variations of the 
fundamental constants. We have derived limits on the 
quadratic interactions of (j) with the photon, electron and 
light quarks from measurements of the primordial "^He 
abundance produced during BBN and atomic dysprosium 
spectroscopy measurements. These limits improve on ex¬ 
isting constraints by up to 15 orders of magnitude. We 
have also derived limits on the previously unconstrained 


linear and quadratic interactions of cj) with the massive 
vector bosons from measurements of the primordial "^He 
abundance. Future laboratory experiments with atomic 
clocks, highly-charged ions, molecules and nuclear clocks 
(see e.g. Refs. [niEi] for overviews of possible systems), 
as well as laser interferometers [59l [60] offer a route for 
exploring new as-yet-unconstrained regions of physical 
parameter space. 
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FIG. 1. (Color online) From top left to right: Limits on the quadratic interactions of (f> with the photon, electron, light quarks 
and massive vector bosons, as functions of the scalar particle mass m,/,. Region below blue line corresponds to constraints derived 
in the present work from consideration of the primordial ^He abundance produced during BBN (the shape of the constraints on 
different sides of the mass ~ 10“^® eV arises from the different dependence of the scalar field amplitude ((>0 on in these 
two limiting regions: (j)o oc => A'x oc for ^ 10“^® eV, whereas (j)o oc => A'x oc for ^ 10“^® 

eV). Region below yelow line corresponds to constraints derived in the present work from consideration of CMB angular power 
spectrum measurements. Region below red line corresponds to constraints derived in the present work using recent atomic 
dysprosium spectroscopy data of Ref. m- Region below black line corresponds to existing constraints from consideration of 
supernova energy loss bounds and fifth-force experimental searches [^. The light quark and massive vector boson interaction 
parameters are defined as (A^)^ = |(A;,)^(A(j)^(md-mu)/[(A(,)^md-(A(i)^m„]| and (A(/)^ = |(AV)^(A'z)^/[(A2)^-1.4(A(v)^ 1I, 
assuming = k'^ and = k'^- 













